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ABSTRACT 

In this paper, we report the detection of a series of radio recombination lines (RRLs) 
in absorption near 26 MHz arising from the largest bound carbon atoms detected 
in space. These atoms, which are more than a million times larger than the ground 
state atoms are undergoing delta transitions (n ~ 1009, An = 4) in the cool tenuous 
medium located in the Perseus arm in front of the supernova remnant, Cassiopeia 
A. Theoretical estimates had shown that atoms which recombined in tenuous media 
are stable up to quantum levels n ~ 1500. Our data indicates that we have detected 
y^Q ' radiation from atoms in states very close to this theoretical limit. We also report high 

, signal-to-noise detections of a, (3 and 7 transitions in carbon atoms arising in the 

■ same clouds. In these data, we find that the increase in line widths with quantum 

number (~oc n 5 ) due to pressure and radiation broadening of lines is much gentler 
than expected from existing models which assume a power law background radiation 
field. This discrepancy had also been noted earlier. The model line widths had been 
^—j I overestimated since the turnover in radiation field of Cassiopeia A at low frequencies 

CZ2 . had been ignored. In this paper, we show that, once the spectral turnover is included 

in the modeling, the slower increase in line width with quantum number is naturally 
explained. 
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1 INTRODUCTION 



Over the forty years after their first detection, RRLs 
have become an im portant probe for interst ellar 
plasm a investigations ijCordoi^^^^^^Sorochenk^^^l 
2002 ). The detection iKonovalenko A. A. fc Sodin L.G.I 



19801: iBlake D.H.. Crutcher R.M. fc Watson W.DJ Il980t 
Koliovatenko^^^^^o^hil7GT^^ir of the low frequency 
RRLs of carbon in absorption towards the strong supernova 
remnant Cassiopeia A (Cas A) in 1980 opened new ways 
to study the tenuous and cool interstellar me dium ( ISM). 
Theo retically, bound atoms with n ~ 1500 iShaver P.A.I 
19751) are expected to exist in tenuous media. Alkali and 
alkaline earth metals like potassium and barium have 
been excited to high quantum levels (up to n ~ 1, 100) 
in the laboratory by using tunable dye lasers and re- 
ducing the Stark effect by minimizing local electric field 
iConnerade. J.-Pl Il998fl . Naturally, the question arose - 
could such large atoms be observed with radio astro- 
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nomical methods? Since atoms in such high quantum 
states trace entirely different ionized regions compared 
to HII regions, they are also important for studying 
physical processes in the ISM. RRL observations till date 
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2001 IPavne H.E.. Anantharamaiah K.R~fc Erickson W.C 
1989t Isorochenko R.L. fc Walmslev C.Mj Il99 ll) cover the 
range from n=166 to n=868 towards Cas A with the 
highest bound state of n = 868 near 20 MHz obse r ved by 
iKonovalenko A.A.. Stepkin S.V. fc Shalunov D.vl 12003) 
using the UTR-2 telescope near Kharkov in Ukraine. At 
such low radio frequencies, terrestrial radio frequency 
interference has to be excised and data collected over many 
days to obtain a high signal to noise detection of any 
feature since the system temperature is fairly high due to 
the non-thermal nature of the background radiation field. In 
addition, the lines broaden due to non-LTE (local thermo- 
dynamical equilibrium) effects and the frequency separation 
between features (1/n 4 ) reduces at low frequencies making 
the observations difficult. 
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In spite of all these difficulties, low frequency RRLs are 
being observed in many Galactic regions. In this paper, we 
report the detection of the largest bound atom in space with 
transitions near n ~ 1009 near 26 MHz towards Cas A. 

While lines at frequencies above 110 
MHz are observed in emission, which im- 
plies the importance of stimulated emi ssion 
jjPa yne H.E.. Anantharamaiah K.R. fc Erickson W.C.I 
n^89F ~the lower frequency features appear in absorption 
against the radiation field of Cas A. Since RRLs have 
been detected at a range of frequencies from the clouds 
in the Perseus arm in front of Cas A, these have been 
extensively used to model and understand the physical 
conditions in the line forming clouds. Since Cas A dom- 
inates the system temperature at all these frequencies, 
it has been easy to use data from widely different tele- 
scopes. It is now widely believed that these low frequency 
RRLs arise in clouds associated with the HI component 
of t he interstellar medium with temperatures of 50 — 75 
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Pavne H.E., Anantharamaiah K.R. & Erickson W.C 
^rochenko R.L. fc Walmslev C.M.I Il99ll) . However, 
most of the models which well explain the variation in the 
line optical depth with quantum numbers, overestimate the 
line widths at low frequencies. 

In this paper, we explain the discrepancy between the 
observations and the existing models used to derive the line 
widths. In section |5] we elaborate on the observations, data 
analysis and present the spectra. In section [3] we discuss 
our results and the width discrepancy and in section ^] we 
present the conclusions of our study. 



2 OBSERVATIONS 

The RRL observations towards Cas A were conducted using 
the biggest decam etre wavelength rad io telescope UTR-2 sit- 
uated in Ukraine jBraude et alll978l> . UTR-2 is a T-shaped 
antenna consisting of 2040 fat dipoles. It has an effective col- 
lecting area of about 100,000 m 2 which results in a signal 
from Cas A which is ~ 10 dB higher than the Galactic back- 
ground at 26 MHz. The telescope has an electronic steering 
capability which can track upto ±70° about the zenith. 

The observations were carried out between 2000 and 
2004 for approximately 100 days. A 4096-channel autocor- 
relometer operating with a bandwidth of about 1.2 MHz 
served as the backend. 

In our case of very high atom quantum states the condi- 
tion of n » An prevails and hence adjacent features were 
considered to be equivalent. Thus, we could fold individual 
transitions and improve measurement sensitivity consider- 
ably. Moreover, such an approach provides more reliable line 
parameters because different transitions were observed si- 
multaneously with the same antenna and ambient medium 
conditions. Fig0 shows the resulting spectrum measured in 
the direction of Cas A in the frequency band around 26 MHz 
with UTR-2. The integration time is 504 h. In the result ob- 
tained from each individual session, four-minute frames were 
cleared from interferences by subtracting from the measured 
autocorrelation function, cosine functions whose parameters 
correspond to hindering signals revealed in the power spec- 
trum. Most often their numbers were between 4-8. Frames 
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Figure 1. Spectrum obtained in the direction of Cas A with 
UTR-2. Relative intensity means ratio Tl/T c . a, (3 and 7 lines 
are visible. 



with wide-band interferences (which were rare) were ex- 
cluded from analysis. Base line correction was carried out 
using splines of the third order. Unblended a and f3 features 
were used to obtain the initial models of the corresponding 
line profiles. After removing of a and (5 lines from the initial 
spectrum and new base line correction we obtained the 7 
lines and were able to build the corresponding model. After 
this, all obtained models (a, /3, and 7 ones) were again sub- 
tracted from the initial spectrum and base line correction 
repeated. Models of the features were built by fitting Voigt 
profiles to folded transitions using Levenberg - Marquardt 
method. The S line model was obtained similarly to the 7 
one. For a and f3 profile fitting the two component model 
with central radial velocities Vi sr (comparatively to the lo- 
cal standard of rest) of -1.6 and -47.2 kms -1 was used. The 
7, and 5 features were fitted with one component profile (at 
-47.2 kms -1 ). The values of central radial velocities as well 
as the value of Doppler width AVd of 15.8 kms -1 were ob- 
tained by two Voigt profiles fitting to the folded a line and 
were fixed when modeling the other features. The procedure 
described above was carried out several times until the fit- 
ted parameters were stabilized. The whole numbers of these 
steps came to six. The Voigt models were used as the ra- 
diative and pressure mechanisms lead to Lorentzian profile 
and the Doppler thermal broadening to a Gaussian. The re- 
sulting line shape is described by the Voigt function which is 
the convolution of the above mentioned curves. Fig|2]reveals 
series of a, f3, 7, and S features. They were obtained as a 
consequent of subtracting the corresponding model profiles 
of RRLs from the original spectrum (see Fig. . 

The individual transitions in Fig.|5|are folded to obtain 
the averaged a, j3, 7 and 8 RRLs of carbon shown in Fig[3] 

The feature near —45 km s _1 is due to gas in the Perseus 
arm. In the a-line spectrum the spectral feature near km/s 
corresponds to the medium in the Orion arm. It is blended 
with the line arising in the Perseus arm and is not the topic 
of consideration in this paper. The effective integration time 
for a transition is 5040 h and radial velocity resolution is 9.2 
fcms -1 . 

The Voigt fitting parameters of lines arising in the 
Perseus and Orion arms are presented in Table 1 (for Orion 
arm features AVd = 4.4 ± 1.1 kms -1 ). 



Radio recombination lines from the largest bound atoms in space 3 

Table 1. Voigt fit to the observed line profiles in Perseus arm 



Galactic arm 


Lines 


T L /T C 


AV L 


AVt 


I 


Perseus 


a(627-636) 


-0.0061 ± 0.00003 


29.0 ±2.9 


36.2 ±3.6 


-28.2 ±2.8 




,9(790-802) 


-0.00088 ± 0.00003 


64.1 ± 10.7 


68.0 ± 11.5 


-8.0 ± 1.5 




7(904-917) 


-0.00032 ± 0.00003 


145.1 ±32.6 


147.0 ±33.1 


-6.3 ± 1.4 




5(994-1009) 


-0.000086 ± 0.00003 


221.5 ±74.6 


222.8 ±75.1 


-2.9 ± 1.0 


Orion 


»(627-636) 


-0.00044 ± 0.00003 


16.2 ±3.5 


17.37 ±3.6 


-0.9 ±0.8 




,9(790-802) 


-0.00014 ± 0.00003 


48.2 ±8.4 


48.7 ±8.5 


-0.8 ±0.8 



Symbols used: Tl/T c = ratio of line temperature to that of background (relative intensity); AVl= Lorentzian line 
width (km s" 1 ); AVt = total profile line width (km s -1 ); I = integral intensity (Hz). 
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Figure 2. The series of a, /3, 7, and <5 radio recombination 
lines observed around 26 MHz in the direction of Cas A. Panel 
(a) shows a series C627.. C636, panel (b) shows series C790 .. 
C802, panel (c) shows 7 series C904 .. C917, and the panel (d) 
shows S series C994 .. C1009. 



3 DISCUSSION 

Previously, the physical models of the medium in the 
direction of Cas A where low frequency carbon spectral 
lines arise have been built based on the analysis of the 
multi frequency data of a transitions because correspond- 
ing features are the most intensive and easy to ob serve 
JPjwn^Lj5^^i<mtfiaraincriah^<^ 

| l994|:lKantharia N.G.. Anantharamaiah K.R. fc Pavne H.E.I 
1998). The model which best fits the data has an electron 
temperature T e =75 K, electron concentration n e = 0.02 
cm -3 and a background radiation field Tb g = 3200 K 
at 100 MHz ( spectral index = —2.6). Dielectroni c-like 
recombination dWalmslev CM. fc Watson W.D.lll982T) and 
the boundary co ndition that the departure coeffi cient 
6 n — » 1 as n — > 00 (Sale m~M. fc Brocklehurst MJI19791) was 
used for the modeling. 

In order to obtain the maximum excitation level of 
an atom which can be studied with the radio astronomical 
methods, we have to analyse the problem of low frequency 
RRLs broadening. This is also important for explanation 
of line parameter behaviour. Pressure and radiative mecha- 
nisms are responsible for the effect along with the Doppler 
thermal broadening. The Lorentzian line width has two com- 
ponents 



AV L = AV P + AVr, 



(1) 



where AVp and AVr are the pressure and the radia- 
tion broadening correspondingly. Towards Cas A radiation 
broadening is caused by influence of the Galactic background 
with temperature Tb g and non-thermal radiation of Cas A 
with WvTcasA (W v is dilution factor) and TcasA is the tem- 
perature of Cas A radiation, i.e 



AVr = AV Rbg + AVncasA- 



(2) 



The theoretic al behaviours of th ese phenomena have already 
been studied dShaver P.AJI197FJ) . The radiation broadening 
is determined by 



AVr 



-1 

£ 

m— 1 



2c/nu > I„B n , m , kms 



(3) 



where I v is intensity of background radiation and B n , m is 
Einstein coefficient of induced emission, c is light veloc- 
ity in km s -1 and v is the observation frequency in H z 
l|Cordon M.A. fc Sorochenko R.LjEool IShaver P. A. II 19 7ft ). 

All models till date have used a power law of index —2.6 
to derive the background temperature. In this case 
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Figure 3. The averaged a, (3, 7 and <5 transitions. Panel (a) 
shows folded a transitions C627.. C636, panel (b) shows folded 
f} transitions C790 .. C802, panel (c) shows folded 7 transitions 
C904 .. C917, and the panel (d) shows folded 6 transitions C994 
.. C1009. 



AVr = 8 * W~ 20 T t n 5 - 8 c/u, kms' 1 



(4) 



where Tt = Tb g + TcasAWu and Tt g is the radi- 
ation temperature of the background emission at 
100 MHz (spectral index = -2.6) iShaver P.A.I 
[ l9_7_5j). When calculating the Lorentzian line w idth 
jKantharia N.G.. Anantharamaiah K.R. fc Pavne H.E.I 
|p99^) ~the total temperature leading to the radiation broad- 
ening was taken as T t = 3200if at 100 MHz (T bg = 800A") 
with spectral index -2.6. This, we believe, is not correct for 
quantum numbers > 650 or so. Such models give a large 
discrepancy between calculated values and experimental 
data. We resolve the problem as follows. Cas A has a 



known turnover point in its spectrum at frequencies less 
than 20 MHz. The flux density drops feraude et allll97Sft 
from 58,000 Jy at 20 MHz to 26,000 Jy at 10 MHz and 
hence instead of increasing with a decrease in frequency, 
the flux density decreases. Earlier models had not taken 
into account this turnover in the spectrum of CasA and 
hence overestimated the line widths at lower frequencies. 
Here we also have taken Tb g = 8007^ when calculating the 
contribution of Galactic background to the total line width. 
In order to determine the contribution of Cas A at the low 
frequencies under study here we have used the experimental 
data for frequency range of 10 - 20 MHz iBraude et alJ 
Il978l) and extrapolated them up to 5 MHz. The dependence 
of flux density and temperature (taking into account the 
dilution factor) can be approximated with good accuracy 
in this range by 



ScasA oc v 



TcasAWv OC V~ 



(•>) 



To evaluate the W v Tc aa A we have used the measurement 
result for n around 630 (y around 26 MHz). Then we calcu- 
lated numerically the values of Ai/RCasA for other frequen- 
cies according to (3). In Fig. 4 this result is graphically pre- 
sented. The dependence of line width against quantum level 
using our algorithm (thick line) and earlier models (thin line) 
along with the experimental data are presented. The param- 
eters are taken as the following: T e =75 K, n e = 0.02 cm' 3 , 
the temperature of the Galactic non-thermal background ra- 
diation T bg = 27,000 K at 26 MHz with spectral index -2.6 
and corresponding dilution factor W v =1 (the drop of Galac- 
tic background is expected to be at v < 4 MHz), and the 
term determining contribution of Cas A emission to radia- 
tive line broadening is estimated as W v Tc as A ~ 12,000 K at 
26 MHz (which is comparable to the influence of Galactic 
background) . ^^^^^^^^^^^ 

The pressure broadening JShaver P.A.Il975l) can be cal- 
culated using 



AV V = 2 * 10" 8 e 



8 -26/T e 



n e n 



2 /Te' 5 c/v,kms . 



(6) 



Correspondence between the curve predicted by our 
arguments and measured parameters is good, whereas the 
older models overestimate the line widths. When we take 
into account the typical conditions of the ISM in the Galaxy, 
the fundamental limit of the highest bound state of a Ry- 
dberg atom is near n= 1,700. These are arrived at by esti- 
mating the quantum level at which the rate of depopulation 
exceeds the orbital frequency of the electr on (Sh aver P. A. 
Il975t iBrocklehurst M. fc Seaton M.j]|l972l) . Our detection 
of the highest bound level at n= 1,005 (1,009) is thus well 
within this range and close to the value quoted earlier 
JShaver P.A]|l975D . This carbon atom is expected to have 
a classical diameter of about 108 fj, ( 0.1 mm). Compare 
this with the Bohr diameter 0.000106 fj.. The carbon atom 
in this tenuous nebula is larger by a factor of ~ 10 6 com- 
pared to the ground state atom. Taking into consideration 
the obtained data and criterion that the distance between 
adjacent lines have to be more than half the line width, we 
can expect that it will be possible to observe the space atoms 
with n=l,100 - 1,200 in the frequency range 15 - 20 MHz. 
Line blending is expected at larger n because of the decreas- 
ing distance between the lines in correspondence with the 
expression: 
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Figure 4. The dependence of line width against quantum num- 
ber. The narrow line corresponds to the best previous model. The 
wide line depicts our model which takes into account the drop of 
Cas A flux density at v < 20 MHz. Points depict experimental 
results. Error bars were calculated using the noise statistics of the 
data. 



Au « (6R c Z 2 An)/n 4 



3v/n 



(7) 



where Rc is the Rydberg constant for carbon and Z is the 
charge of an ion. 



4 CONCLUSIONS 

In this paper, we have reported the detection of the largest 
bound atom in space, which has recombined to quantum 
levels of ~ 1009. The final spectrum presented here shows 
the presence of a, /3, 7 and 8 lines. This spectrum near 
26 MHz was obtained by integrating data towards Cas A 
obtained using UTR-2 telescope in Ukraine. 

At such low frequencies, the line widths are expected 
to increase steeply (~oc n 5 ) due to radiation and pressure 
broadening which strongly dominate over thermal Doppler 
broadening. However, we find that the observed increase is 
much slower than predicted. We believe that this slow in- 
crease is due to the decrease in the background radiation 
field at low frequencies due to the spectral turnover in the 
non-thermal spectra of Cas A. After incorporating this in 
the model, we find an excellent fit to the line widths of our 
observed data. Thus, we suggest that this decrease in the 
field should be incorporated in any model which explains 
the RRL data towards Cas A. 

Presently investigations of low frequency RRLs with 
UTR-2 are continuing for many galactic objects and data 
processing is progressing. It is clear that RRLs studies will 



be important part of the scientific programs of future giant 
low frequency instruments such as LOFAR. 
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